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Abstract 

We present a fast hadronization model for the constituent quark plasma (CQP) 
produced in relativistic heavy ion collisions at SPS. The model is based on rate 
equations and on a confining equation of state inspired by the string phenomenol- 
ogy. We display the time evolution of the relevant physical quantities during the 
hadronization process and the final hadron multiplicities. The results indicate that 
the hadronization of CQP is fast. 

1 INTRODUCTION 

The research of ultra-relativistic heavy ion collisions is from its beginning motivated by 
the search for quark matter: a bulk of deconfined color charges. The prime attempts to 
describe the hadronization assumed, that this quark matter consists of non-interacting 
massless quarks and gluons [1], 0, [J. In place of these early image of free quarks and 
gluons (plasma) - which was based on thermodynamical studies of the pure nonabelian 
gauge theory on the lattice - gradually the picture of a quark matter emerges, containing 
effective propagators and interaction vertices jjjj, j|, ||. At the characteristic energy scale 
of CERN SPS heavy ion experiments (y/s/2 w 10 GeV/nucleon) the dressed gluons are 
heavier than the quarks (M g « 600 - 800 MeV, M q « 150 - 300 MeV) [g, and both 
values are bigger than the temperature (T « 160 MeV). Therefore the number of quarks 
overweight that of the gluons and can be treated non- relativist ically. On the ground 
of these theoretical indications we expect that in the CERN heavy ion experiments not 
an ideal quark gluon plasma, but a constituent quark plasma (CQP ), which contains 
antiquarks too, is formed in some intermediate state of the reaction (see the ALCOR 
model, Ref. [7|, §). 

For the description of heavy ion reaction products at SPS a purely hadronic interpre- 
tation of bulk experimental results was also suggested |J. This alternative interpretation, 
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however, has problems in explaining the process of the creation of new hadrons, especially 
the short time scale. Namely, the hadronic processes have characteristic times of several 
tens of fm/c |3], TO] , while the typical heavy ion reaction time at SPS energy is about 1—2 



fm/c [jy]]. Processes on the quark level on the other hand have equilibration times of the 
order of 0.1 - 2 fm/c |g0. 

The CQP picture has some advantages in comparison with the quark-gluon plasma 
scenario too: In the hadronization of a quark-gluon plasma the principles of color con- 
finement and entropy generation can be both satisfied only with the assumption of an 



extreme growth of the reaction volume |13[ . Further, a slow first order transition through 
near-equilibrium states would need too long time (up to 50 fm/c), due to the re-heating 
TJ. 



In this paper we assume that a massive quark matter, the CQP, is formed in heavy 
ion reactions at SPS energy, which is in thermal but not necessarily in chemical equilib- 
rium. The evolution of this CQP is followed through a set of coupled time dependent 
differential equations for the temperature and chemical composition. The construction of 
the equations requires relations between the different equilibrium particle numbers at an 
actual temperature as input. These input values are determined by the equation of state. 

We mention in this place, that the ALCOR model ]7| §], the hadronization problem 
is treated with an algebraic approximation, in contrast to the method of the present paper, 
where we follow the complete time evolution of the system. The inclusion of confining 
equation of state in the present model is very important difference to the ALCOR model, 
and to our best knowledge this type of phenomenological equation of state with confining 
character has not been discussed in the literature until now. 

The paper is structured as follows: in Section |[ we discuss the equation of state 
of the mixture of hadrons and interacting quarks. In Section [3|. we describe in details 
the dynamics of the hadronization. In Section |j. we discuss our numerical results. The 
conclusion is drawn in Section^. In Appendix A. we collected the relevant hydrodynamical 
and thermodynamical and statistical expressions used in the paper. 



2 EQUATION OF STATE OF THE MIXTURE OF 
HADRONS AND INTERACTING QUARKS 

We assume that at the beginning of the hadronization the matter consists of massive 
quarks and anti-quarks. In the time evolution of the system quarks and anti-quarks 
form diquarks, anti-diquarks, mesons, baryons and anti-baryons. We assume that the 
mixture of all of these particles are in thermal equilibrium which can be characterized 
by a temperature. For the representation of the interaction of the colored particles we 
introduce an extra term into the free energy, which is inspired by the string picture. 

We note that there is an important difference to the canonical approach to color 
confinement transition: in an (ideal) mixture of quarks and hadrons the occupied volume, 
V, is the same for both components, V q = Vu — V, the pressure contributions p q and ph 
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are additive. On the other hand in the application of the Gibbs criteria of a phase co- 
existence the volumes V q and Vh are additive V q + Vh = constant, and the partial pressures 
are equal, p q = Ph = P, in phase-equilibrium. In our physical picture of hadronization 
there is no phase coexistence and the Gibbs criteria do not apply. Colored particles 
and color-neutral clusters, pre-hadrons are distributed in a common reaction volume and 
chemical reactions convert eventually the quark matter into a pure hadronic matter. 

During this hadronization process the reaction volume expands and cools. Due to 
the change of the multi-particle composition this expansion is not adiabatic: some heat 
can be produced (or consumed) by quarkochemical processes. The expansion law of an 
ideal mixture follows from eq. fl44f) and eq. (|59|) in Appendix A: 

E + \tY, Ni + §T£ K + £ N< = 0. (1) 

i i i t 

Due to the foregoing hadronization the number of particles decreases, 

E^<°> (2) 

i 

therefore this process re-heats the system. Cooling effects are due to the expansion 
(V /V = <9 M w M > 0) and rest mass creation (J2i m iNi > 0). The latter is possible with sat- 
isfying (0) only if hadron masses are larger than the sum of their constituents. 

At the present level, however, the physical process of color confinement by the 
hadronization is not yet taken into account. In particular due to the dilution during the 
expansion it happens that not all quarks or diquarks find a partner to hadronize. In order 
to avoid this effect we supplement the model by the following confinement principle: all 
particles carrying color charge (quarks, diquarks, anti-quarks and anti-diquarks) will be 
penalized by a free energy contribution stemming from strings. The number of strings is 
proportional to a weighted sum of the number of color charges, 

Q = E Qi N i, (3) 

i 

Here it is g» = for hadrons, g« = 1/2 for quarks and anti-quarks, and g« = 3/4 for 
diquarks and anti-diquarks. The higher effective charge of diquarks reflects a possibly 
higher number of in-medium partners, to which a string is stretched. 

The average length L of a string depends on the density of colored objects, as 
L = n^ 1 / 3 , where 

n c = Y,Ni/V : (4) 

where the summation i 6 c excludes color neutral particles (hadrons). So the free energy 
of the ideal quark matter - hadron matter mixture, F = F^ + AF, is supplemented by 
the following contribution of strings: 

AF = a s n: 1/3 Q, (5) 
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with the effective string tension a s ~ 1.0 GeV/fm. We note here that if all colored 
particles carried the same color charge, the effective string tension would be equal for all, 
and the interaction free energy would become 

AF = a s , eS n 2 J 3 V. (6) 

In order to achieve baryon production and diquark elimination properly we shall, however, 
need to use the more complicated ansatz eq.(|5]). 

This additional free energy comprises the non-ideality of the equation of state we 
use. Since this addition is proportional to the volume V and the rest depends on densities 
only, it satisfies thermodynamical consistency requirements |15| due to its construction. 

While there is no new contribution to the entropy, 

S = S id , (7) 

the pressure, the energy and the chemical potentials of colored (qi ^ 0) particles receive 
important modifications: 



P = Pid~ ^v s n c 



i -1/3Q 

3 c V 
E = E iA + a s nZ 1/z Q, 

Hi = A*i,id + (r s n~ 1/s (qi - -q), (8) 

with q = Q/(Vn c ). Hadronic chemical potentials have no modifications at all. 

This non-ideal completion of the equation of state influences both the expansion 
and cooling and the changes of particle composition. Since the entropy is not changed by 
the introduction of the interaction term (see. eq. (|7|)) therefore 

TdS = TdS iA . (9) 

From this equation, using eq. (|46|) we obtain 

dE + pdV-Y, HidNi = {dE + P dV) id - J2 Vi,iddNi . (10) 



Applying eq. (fH|), the non-ideal cooling law becomes 

(dE + P dV) id + £ (/ii - m M ) dN i = 0. (11) 

i 

Accordingly eq.([I|) is supplemented by a new generic term due to the non-ideal equation 
of state, 

T 2 V EiN t 2 Z i (m i /T)N i 

T 3V EiN t 3 

2 Ei(Hi/T - fi iM /T)Ni 

3 EiNi 1 } 
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The additional term has the form 

E ^p^N, = J> - ^)N t . (13) 

This term is negative if color charges are eliminated from the mixture. Therefore color 
confinement, causing an extra suppression of equilibrium numbers of quarks and alike 
particles, re-heats the expanding fireball as well as the "normal" chemistry of the ideal 
quark - hadron mixture. The more one suppresses color charges with respect to an ideal 
mixture the more re-heating occurs during hadronization. This, of course, works against 
the hadronization process. The only physical effect besides a fast expansion - which has, 
however, kinematical limits stemming from scaling relativistic expansion - that can cool 
the mixture sufficiently is rest-mass production. The color charge eliminating hadroniza- 
tion therefore must be accompanied by the production of heavy hadron resonances. (Rest 
mass production due to quark pair or gluon creation would not reduce color.) Since in our 
equation of state we have an explicit interaction energy between the quarks, our effective 
quark masses should be less than that given in Ref . || . We shall use the following values: 
(m u = [m 2 + m^] 1 / 2 , = [m^ + m 2 ^] 1 / 2 , m s = [m 2 + m 2 ^] 1 / 2 , with thermal mass 
m th = 0.15 GeV and m u0 = m d0 « 0, m s0 = 0.15 GeV. The clusters have a mass accord- 
ing to the average mass extra to the summed valence quark masses of the two lowest lying 
hadron multiplets: the pseudoscalar and vector meson nonets and the baryon octet and 
decuplet, respectively. 



3 DYNAMICS OF HADRONIZATION 

3.1 Initial state 

The initial energy density — distributed along the beam direction between — To sinh^o an d 
tq sinh 7]o — can be related to the center of mass bombarding energy y^i in the experiment, 



£ ° nR 2 r 2sh m (14) 
On the other hand the initial invariant volume dual to dr at constant r = r is given by 

V = nR 2 r 2 Vo . (15) 

The initial internal energy (i.e. the energy without the collective flow of a fluid cell) at 
t = t is therefore less than y/s for finite 770: 

£0 = £ Vq = -^v^- (16) 

At the CERN SPS experiment 770 ~ 1.75 (due to some stopping), R rs 7 fm, r ~ 0.8 
fm and we obtain Vq ~ 431 fm 3 and Eq w 2.13 TeV. Compared to the total energy 
of about yfs = 3.4 TeV (carried by about 390 participating nucleons in a central Pb- 
Pb collision) approximately two third of the energy is invested into rest mass of newly 
produced particles and thermal motion and and one third into the flow. 
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Comparing this with an alternative expression for the thermal energy of an ideal 
massive quark matter, 

i 

one can estimate the initial temperature at the beginning of hadronization. Using our 
standard values for the incoming quark numbers []16 N u (0) = 544, AT d (0) = 626, further, 
assuming that 400 uu, 400 dd and, with f s = 0.21, 168 ss quark anti-quark pairs are 
created in one central collision, we get from the above equation T = 0.18 GeV. We use 
these numbers for the newly produced quark pairs in order to arrive the experimentally 
measured hadron and strange particle numbers. 



3.2 Hadronization processes 

Chemical equilibrium is not supposed initially, rather a definite over-saturation of quarks 
in the reaction volume. The initially missing color-neutral hadron states - mesons and 
baryons - are formed due to quark fusion processes in a non-relativistic Coulomb potential. 
The rates for different flavor compositions differ mainly due to the different reduced masses 
of quark anti-quark or quark diquark pairs. First of all this influences the Bohr radius 
in the Coulomb potential J7|. Of course, the presence of a medium - which establishes 
the necessary momentum balance after the fusion - also influences the hadronization rate. 
The cross section for such a 2 —>■ 1 process in medium is 



p\ 3 16M 2 V7ra 
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a = - v - ^ 18 

with a = l/(am) Bohr radius of the Is state in the Coulomb potential and p is the Debye 
screening length [|7| . Here p is the relative momentum of the hadronizing precursors, m 
is their reduced mass and M is the total mass. 

The coupling constant we are using is a function of the relative momentum according 
to the formula 

-^r = — + /Wog-^- (19) 
or(p) a A 2 QC b 

as long as a(p) < atQ. At smaller momenta it levels off at «o- Here 

^ cd = h (t " 2Nf ) m 

is the beta-parameter of the one loop beta-function occurring in the charge renormaliza- 
tion process in perturbative QCD. The function (|T^) is a good phenomenological approx- 
imation and has the correct infrared and ultraviolet asymptotics. The value cuo ~ 1.4 is 
taken in order to fit the gluon condensate strength in vacuum. 

The medium effect is parameterized in all cross sections by p ~ 0.2 — 0.3 fm in 
a flavor independent way. If the hadronization process is fast enough, the prehadron 
numbers (i.e. hadron numbers before resonance decay) become proportional to the two- 
body reaction rates. Since quark fusion does not change the number of valence quarks 
during the hadronization the fast hadronization limit leads to the ALCOR model . 
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The relative momenta are taken from a random Gaussian distribution, 

dV(p) oc e~ p2 / 2mT d?p, (21) 

at temperature T and reduced mass m. This method allows us both to simulate thermally 
averaged hadronization rates, 

J2 =<"—>, (22) 
m 

or to follow the event by event variation of the quark - hadron composition. For a fast 
estimate of thermally averaged rates about 3 — 5 points in the relative phase space suffice 
for each reaction at each time instant. 

Although the above corrections to the equation of state help to reduce the equilib- 
rium ratio of free quarks to those confined in hadrons, total color confinement would only 
occur after a long time when the mixture is cool and dilute. This is obviously not the 
case in relativistic heavy ion collisions. Therefore we also apply a dynamic confinement 
mechanism in our model: the medium screening length p occurring in the hadronization 
cross section will be related to strings pulled by color charges trying to leave the reaction 
zone. This way the screening length p is increased as the color density decreases: we keep, 
however, the product p 3 N c /V constant, 

p(t)=p(0)*(N c (0)/NM) 1/3 - (23) 

Here N c stands for the number of colored objects (each in triplet or anti-triplet represen- 
tation), 

N c = N q + Nq + N d + Njy. (24) 

As we shall discuss in the next chapter both color confinement mechanisms are necessary in 
order to achieve a pure hadronic composition in a short time by trans-chemical processes, 
while satisfying the requirements of the entropy growth and energy conservation. 



3.3 Reaction network 

What remains to specify the model is the system of rate equations describing the trans- 
formation of quark matter into hadronic matter. We consider Np = 3 light quark flavors 
u, d and s. 

There are Np(Np + l)/2 = 6 possible diquark flavors and the same number of anti- 
diquark flavors. The number of quark anti-quark flavor combinations is Np = 9 while 
that of quark or anti-quark triplet combinations is N F (N F + 1)(N F + 2)/6 = 10. In the 
hadronizing quark matter we deal with altogether 2*3 + 2*6 + 9 + 2* 10 = 47 sorts of 
particles. 

Let us generally denote quarks by Q, diquarks by D, mesons by M and baryons by 
B. The possible fusion reactions are: 

Q + Q — > D, 
Q + Q — > D, 
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Q + Q — > M, 
Q + D — > 5, 

Q + £ — ► 5. (25) 

Accordingly the number of reaction channels, 2* 10 + 9 + 2*6 = 41 is 2iVV = 6 less than 
the number of particle sorts. Therefore in chemical equilibrium, when the above reactions 
are balanced by the respective decays, the number of all particles are under-determined. 

This is fortunately not a problem, because the hadronization reactions under con- 
sideration do not create or annihilate elementary quark flavors, they just produce new 
combinations. There are therefore exactly 2Np = 6 quark and anti-quark flavor numbers, 
which are conserved by these reactions. They are determined by the initial state. 



Our model is completed by the system 
reaction of type 

i+j — 

we account for the changes 

dNi = dNj = -Adt, 
cumulatively in each reaction. Here 



of rate equations. Considering a general 
-> k 

dN k = +Adt, (26) 



A = Ri^tNiNj (l - (27) 

with a thermally averaged rate Rij^k (cf. eq.(^)). The changes stemming from different 
reactions accumulate to a total change of each particle sort in a time-step dt. 

Inspecting the expression ( p7|) it is transparent that an extra increase of a chemical 
potential suppresses the production or enhances the decay of the corresponding particle. 
This behavior ensures that 

S = -T,%Ni>0, (28) 

i 

for any complicated network of reactions. The approach to chemical equilibrium always 
produces entropy. 



3.4 Chemical equilibrium 



Chemical equilibrium is defined by the requirement that all chemical rates vanish (cf. 
eq.^7]). It leads to relations like 

/4 q + »7 = ( 29 ) 
for each reaction channel. The correspondence between equilibrium chemical potentials 
and equilibrium number densities is, however, in the general case not as simple as for a 
mixture of ideal gases. From eq.(|J) we obtain an implicit equation for the equilibrium 
densities of colored particles, 



n 



cq 



nf 1 exp 



/V 1 - bi(nl q ) - nii 



(30) 
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where n^ q is the color charge density in equilibrium, (see Eq.(H)), nf 1 = N\ h /V from 
Eq.© and 

b l (n c ) = a s n- 1 l*(q l - l -q} . (31) 

We call the attention to the fact that no solution of Eq. (|3H|) exists below a critical tem- 
perature. 



Since the number of reactions is less than the number of particle sorts, the equilib- 
rium state is not fully determined by these conditions alone. The missing information is 
contained in the value of conserved numbers. 

Applying eq.(|8|), the non-equilibrium chemical potentials, and hence the essential 
factors, e~ Mi//T , in the chemical rates (eq.(^7])) can be expressed as 

e = e *i l 1 . e T . -±-. (32) 

In the combinations appearing in the detailed balance factor of the rate equations using 
eq . (|29|) we obtain 

'-^ = 1 -t#~ (33) 

with 

A^ i = 6 i (iV)-6 i (iV^). (34) 

The corrections 6»(iV) in the non-equilibrium chemical potentials may in general depend 
on the number densities of several other components on the mixture. 

At this point we note that the extra e~ Afl ^ T factors occur for non-ideal equations of 
state where the correction to the free energy density is a nonlinear function of the number 
densities. 



3.5 Detailed rate equations 

With the purpose of better understanding in mind we enumerate certain types of the 



rate equations given in the above compressed form in eqs.( |26| , p7|) . Let {u, d, s, . . .} be an 
ordered set of flavor indices. We denote quarks by flavor indices i, j or k, the anti-quarks 
by their overlined versions, diquarks consisting of a flavor pair satisfying i < j by [ij], a 
meson made from a flavor % quark and a flavor j anti-quark by [ij], and finally a baryon 
made from quark flavors i,j and k satisfying i < j < k by [ijk]. In fusion reactions we 
create primarily mesons, (anti) diquarks and (anti)baryons. 

The rate equations for these reactions are straightforward to encode on a computer 
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but lead to somewhat clumsy equations writing down the actual rate equations: 

( N eq N- q 

j \ [iJ] 

( N eq ]\[ eq 

£ R»m-+m NiN w - -L-Miv^e^^-Aw-^j/T I (35) 



N eq 

[jk] V ly [ijk] 



for quarks, 



V M / i<j 

(N eq N eq 
N M N k - J^ Nme ^m-*m-<**rr ] ( 36 ) 
iy [ijk\ 



for diquarks, 

A AT - / N eq N- q 

yJ ir = - -T^f- %^"'«- a "-^ /T ) (37) 

for mesons and finally 

V i*j fc i / j<j<fe 

/ N eq N? q i N 

+ i2,- +[ * My *] W*] - _JJ7V [ .. fe]e (A, fo , 1 -A^-A, [ifc] )/T 

V fo' fc ] / i<j<fc 

/ 7\r e 9 ]y eq \ 

+ Rk+m+m n * n m - ^^AW^* 1- ^ - ^ 1 ^ (38) 

V 7V fo' fe ] / i< j<k 

for baryons. We obtain similar equations for anti-quarks, anti-diquarks and anti-baryons. 



We would like to point out that the novel appearance of the above rate equations 
is a consequence of the equation of state we are using. This form goes beyond the usual 
expression multilinear in iVj/iV i eq , and reflects the fact that the hadronization reactions are 
strongly influenced by the presence of colored particles not taking part in the particular 
two-body process directly. 
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Especially the rate equation for baryon production is complicated. The conditions 
indicated as lower indices at the bracket expressions or in the summation signs rule out 
multiple counting of baryons which contain equal flavor quarks. For example the reaction 
u + [uu] — > [uuu] makes only one [uuu] baryon (pre-A ++ ) and not three. On the other 
hand [uds] can be constructed from three different quark - diquark fusion process 

u + [ds] — ► [uds] , 
d + [us] — > [uds] , 

s+[ud] — ► [uds]. (39) 

The only easy and transparent notation we found for these reactions is the compact 
eq. (|26|) . Since we calculate the evolution of the number of a given type of quark group (e.g. 
uud), at the final time we distribute this number between the corresponding multiplets 
- the lowest lying pseudoscalar and vector nonets for the mesons and lowest lying octet 
and decuplet for baryons (antibaryons) — according to the spin degeneracy. (E.g. the 
number of uud quark group is distributed between the p + and A + in the ratio two to four 
while uuu populates only the A ++ resonance. ) 



3.6 Hadronic decays 

The set of rate equations describes the time evolution of the number of all involved 
particles. In order to get the final hadron numbers we integrate these equations until the 
number of colored particles becomes negligible. This way we obtain a number of hadronic 
resonances (in the present version the vector meson nonet and baryon decuplet). Finally 
hadronic decays are taken into account with the dominant branching ratios obtained from 
Particle Data Table jD]]. We assume that secondary hadron-hadron interactions have a 



negligible effect on the finally observed hadronic composition. 

The time evolution of the entropy and temperature is obtained by simultaneous 
integration of eqs. (|58|) and ([[]). The structure of the rate equations containing both 
creation and decay terms ensures that the entropy never decreases during the transition 
2jfl. On the other hand the energy conservation is established by the effective cooling 



law eq.(|]). The numerical study of this complex system of equations is presented in the 
next section. 



4 Numerical results and discussion 

In this section we present the results of the numerical solution of the set of the 41 + 
2 coupled differential equations. In the Figures the calculated time evolution of the 
temperature, entropy, internal energy and the different particle and antiparticle numbers 
are presented for a 158 GeV / nucleon Pb + Pb central collision. For the numerical 
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solution the initial condition and the values of the parameters has to be specified. For the 
parameters describing the initial state we used those given in subsection |3.1| , while for the 
parameters determining the dynamics of the hadronization we used the following values: 
p = 0.2/m, a = 1.4. Both the quarks and diquarks vanished only in the case, if we used 
the value qa — 1.5 for the effective diquark color string tension. The time evolution of 
different quantities calculated with this initial condition and parameters are shown in the 
Figures. 

From Fig. la we can see, that at the beginning of the hadronization there is a rapid 
decrease in the temperature due to the rest mass formation of the hadrons. Shortly after 
that, the reheating starts as an effect of color confinement (see eqs.(28) and (32)). The 
eq.(28) shows, that by removing two colored object to produce a colorless hadron, the 
associated string energy is also removed, and it has to appear in the thermal energy. 
Since this energy is inversely proportional to the one third root of the color density, this 
effect is stronger at smaller color density. Finally, as the hadronization is completed, 
the expansion bleeds to the cooling of the system. Fig. lb shows, that the total entropy 
is monotonically increasing during the hadronization, while the contribution of colored 
particles (CQP) is gradually eliminated. In Fig. lc one can observe an interesting pattern 
in the time evolution of the pressure. The partial pressure of the interacting CQP rapidly 
decreases as the number of quarks decrease. As the color density drops, this pressure 
becomes even negative. The increasing hadron partial pressure, however, overcompensates 
this negative value. At later times with the expansion of the ideal hadron gas the total 
pressure decreases. Due to this interplay, the total pressure falls steeply in the early 
period, it stays for a short while almost constant in the re-heating period, and finally 
drops slowly at the late expansion stage. Assuming an initial condition leading to slower 
hadronization, this balance is not so effective and in the time evolution of pressure a local 
minimum appears. The partial and total internal energy evolution, displayed in Fig. Id, 
shows, that the hadronization is completed at 2 fm/c after the beginning of the process. 
The decrease of the internal energy is compensated by the work of pressure while making 



the flow (cf. eq.([44])), the total energy density drops due to the volume expansion like 



e oc 1/t at very late times only, when p = 0. 

Fig. 2 shows the time evolution of different colored particles. The quark and anti- 
quarks numbers are monotonically decreasing functions of time. The diquarks are pro- 
duced from the quarks, and then they contribute rapidly to the formation of baryons. 
That keeps their number always on a low level. In Fig. 3 the evolution of hadron numbers 
is shown. Here we did not take into account the decay of the resonances. Actually these 
time evolutions are the evolution of the constituent quark clusters, i.e. prime hadronic 
resonances, from which the experimentally observed particles emerge. It is clear from this 
Figures, that the mesons are formed faster than the baryons. This is understandable, 
if we consider, that the mesons are produced in a one step process, while the baryons 
are formed in a two step process: first the formation of diquarks and in a diquark quark 
reaction the formation of the baryon. In the ALCOR model the ratio of hadronic species 
are determined by the ratio of steepnesses of these curves. Since these curves do not cross 
each other, one can understand, why the algebraic ALCOR approach to the solution of 
rate equations is a good approximation. The difference, that at the very beginning the 
meson curves increase linearly with the time and the baryon curves start with a quadratic 
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form, was taken into account in the ALCOR model by one single common factor, the so 
called baryon suppression factor. 



In our model we calculate the number of all produced hadrons. In Table 1 the 
hadron numbers obtained with the Transchemistry model and those obtained with the 
ALCOR [16] and the RQMD ]18| models are shown together with the few published 



experimental data. The present results of the ALCOR model are slightly different from 
the previously published ones |Tj|, because here we considered the feeding of A (and 
A ) particle multiplicities from S°, H - , E°, Q~ (and S°, S + , H°, Q + respectively), as it 
was obtained in the experiments. The fit parameters to this new case are: N qq = 391, 
N S g = 172 and a e // = 0.97 comparing to the earlier N qq = 398, N S g = 175 and a e ff = 1.03. 

Table 2 shows a comparison for the multi-strange baryon ratios. While in many 
cases there are intriguing agreements, in some other cases there are some discrepancies. 
The reason for this may originate from two sources: i) the experimental data referred 
here are the production ratios in the overlap window of the detector acceptances. Thus, 
if the momentum distribution of the two particle sorts is not the same, then these ratios 
are not equal to the total number ratios, ii) these calculated values are more sensitive to 
the simplifying assumption, that the hadronization happens into the lowest lying baryon 
octet and decuplet and the two lowest lying meson nonet. We have committed this 
simplification, because the inclusion of other hadron multiplets would have multiplied the 
workload in the computer simulation. In our model we don't have such parameters with 
which the different particle ratios could be manipulated independently. The degree of 
the overall harmony between our results and the recent experimental data seems to be 
promising and thus it could initiate the search for further signatures of the CQP state. 



5 Conclusion 

In this paper we presented a new model for the hadronization of a constituent quark 
plasma (CQP) based on rate equations in a quark matter - hadron matter mixture. The 
color confinement was taken into account by using consistently a plausible equation of 
state motivated by the string model. The role of the different physical processes entering 
into the hadronization was discussed. Our results presented in the Figures clearly show 
a very fast hadronization. This is mainly a consequence of the large hadronization rates, 
Rij^k an d of the fact that in the initial state the system is very far from equilibrium. 
Observing the shape of the time evolution of different hadron multiplicities, it became 
understandable, why the simpler algebraic approximation, applied in the ALCOR model, 
works so well. The comparison with the existing experimental data indicate, that it is 
possible, that in the PbPb collision at SPS a piece of matter is formed, inside which the 
massive quarks and anti-quarks interact with a string like mean field. 

Finally we emphasize, that this type of phenomenological investigations are neces- 
sary, as long as the hadronization of the quark matter as a non-equilibrium, non-static, 
non-perturbative process, cannot be described with the methods of the elementary field 
theory. 
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Figure 1. The time evolution of the temperature T (a), the entropy S 
(b), the pressure p (c) and the energy density e (d) of the system together 
with the partial contributions of quarks and hadrons. 
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Figure 2. The time evolution of the number of colored particles. The 
line styles of different flavor compositions are indicated in the respective 
figures. 
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Figure 3. The time evolution of the number of color neutral clusters 
consisting of quarks and/or anti-quarks indicated on the figures. 
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Table 1: Total hadron multiplicities for Pb+Pb collision at 158 GeV/nucleon bombarding 
energy. The displayed experimental results are from the NA49 Collaboration: estimated 
result a is from JHJ; b is from [20]; c is from [^T|; * is estimated from {K~} and {K®}. 
Theoretical results are from the Transchemistry, ALCOR and RQMD |l8j ("ropes + no 
rescattering" version) model. Here it is A — like = A + S° + H~ 
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Pb+Pb 
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ALCOR 


A u /A° 


n i ?8 + n m 2 


Oil? 


117 




266 ± 028 


0.243 


0.227 




0.46 ±0.15 


0.529 


0.564 


E~/A° 


0.093 ±0.007 


0.078 


0.071 




0.195 ±0.023 


0.170 


0.138 


Q-/E- 


0.195 ±0.028 


0.095 


0.125 



Table 2: Strange baryon and anti-baryon ratios measured by WA97 Collaboration p3"| 
and obtained from Transchemistry and ALCOR model for Pb ± Pb collision at 158 
GeV/nucl bombarding energy. The experimental data are the production ratios in mid 
rapidity at pt > GeV. 
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APPENDIX A: 

THE RELEVANT HYDRO DYNAMICAL AND 
THERMODYNAMIC AL EXPRESSIONS 

The familiar energy-momentum conservation 

d v T^ v = 0, (40) 

with 

r = (£+p)!iV-pf, (41) 

being the energy-momentum tensor of a perfect fluid with flow four-velocity u M , local 
energy density e and pressure p, has two interesting projections: one is parallel to the 
four-velocity, 

u^duT^ = d u (eu u ) + p{d v u v ) = 0, (42) 
and the other is orthogonal to this, 

(u^(u v d v )-d^)p + (e+p)(u v d v )u^ = 0. (43) 
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The former - expressing local energy conservation - can be casted into the suggestive form 

dE + pdV = (44) 

using an infinitesimal volume element, dV = (d^u^Vdr, the total internal energy, E = 
eV, and the parametric derivation with respect to the time- variable in the co-moving fluid 
cell, 

4- = u v d v . (45) 

GST 

Starting from the first principal theorem of thermodynamics 

dE = TdS - pdV + Vi dN h (46) 

and using eq.([44]) we arrive at the following change of rate of the total entropy: 

S = "Ef^ (47) 

i 

where the 'dot' stands for the time derivative (S = dS/dt). The chemical potentials 
can generally be derived from the equation of state. 

We consider a flow pattern given by the Bj or ken-flow 

Ufl = (t^0,0). (48) 



Here 



T T 



2 j.2 2 

T = t — Z , 



(49) 

and the four-velocity is normalized to = 1 and its four-divergence is 

dX = ~ ■ (50) 
r 

The energy density e and the pressure p depend only on r and the flow is stationary. 
Corrections stemming from a mild radial flow also can be worked out, but we do not use 
these corrections here. The effective dimensionality of the flow is 1 in this case. 

Since we deal with a system of rate equations for particle species numbers iVj occu- 
pying a common reaction volume V at temperature T, the generating thermodynamical 
functional is the free energy, F(V,T, A^). The chemical potentials are given by 

OF 

" = W' (51) 

the pressure by 

dF 

(52) 
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the entropy by 

dF 

s = -gr> (53) 

and finally the internal energy is 

E = F + TS. (54) 

As a starting point we consider a mixture of ideal gases of massive quarks, diquarks, 
mesons and baryons, and their respective anti-particles. The corresponding free energy is 

F id = ]T TN t ( In - 1 J + m ^ ( 55 ) 

with Maxwell-Boltzmann statistics for the non-relativistic massive matter: 

Af = Vdif-0T 3 e- p2/2m < T (56) 

Here di = (2si + 1)q are spin and color degeneracy factors. The chemical potentials in 
an ideal gas mixture are 

/^id = Tin j-^ + mi . (57) 



The total entropy is given by 



Sid = £A^-m^). (5f 



The energy and pressure of such an ideal, non-relativistic mixture is given by 



£id = Y.[ m i + \ T ) N i 



Pid = (59) 
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